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Abstract 

The Maxwell equations are formulated in a generally covariant and metric- free way in 1+3 
and subsequently in 4 dimensions. For this purpose, we use the excitations T>, Ti and the field 
strengths E,B. A local and linear constitutive law between excitations and field strengths is 
assumed, with a constitutive tensor x^ kl = —y} tkl = —^ lk of 36 components. The properties 
of this tensor are discussed. In particular, we address the validity of the Post constraint, 
a subject that is very much under discussion. In this connection, the Tellegen gyrator, the 
axion field, and the "perfect electromagnetic conductor" of Lindell & Sihvola are compared 
with each other. 
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1 Introduction 

Electrical engineers as well as physicists are familiar with the Maxwell equations 

dV = p, dH-V = j, (1) 
dB = 0, dE + B = 0, (2) 

here expressed [1, 4, 5, 19, 39] in exterior differential forms in terms of the electric and magnetic 
excitations V, 7i and the electric and magnetic field strength E, B, respectively. With d we denote 
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the 3-dimensional exterior derivative and with a dot the time derivative. The sources of the in- 
homogeneous equations are the electric charge density p and the electric current density j. These 
equations are expressed in a generally covariant form, that is, they are valid in arbitrary curvilinear 
coordinates. Moreover, they are metric-free and thus the concept of orthogonality is not used at 
all. 

As shown in the literature, see [40, 34, 5], the Maxwell equations can alternatively be written 
in terms of components as 

d a V a = p, e abc d b n c -V a = f, (3) 
d a B a = , e abc d b E c + B a = , (4) 

with a, b, ■ ■ ■ = 1,2,3 and e abc = ±1,0 as the totally antisymmetric Levi-Civita symbol. Here 
D a , B a ,j a are vector densities, E a , H a covectors, and p is as scalar density. Also in this formulation 
the general covariance and the condition of being free of the metric are manifest. Note, since no 
metric is used, we must not raise and lower indices with the help of the metric at this stage. 

The Maxwell equations need to be supplemented by constitutive relations. And this for the 
description (i) of empty spacetime (the vacuum) and (ii) of material media. We will discuss here 
the simplest case, namely local and linear media. Local in space as well as in time, that is, the 
excitations at every point in space and time depend only on the field strength at the very same 
point. 

A general local and linear constitutive relation reads (in the conventions of [5]), 

H a = -£ b a E b + <B ba B b , (5) 

V a = -%ba Eb + :Db a B b (g) 

The generalized permittivity matrix 2t ba and the generalized impermeability matrix (the reciprocal 
of the permeability) 23& a are well known. Less known, but well established in the electrodynam- 
ics of continuous media, see [31], are the magneto-electric matrices £* a and D b a , which describe 
magnetic-electric cross effects (like the Faraday effect or optical activity). Altogether, these matri- 
ces encompass 4 x 9 = 36 independent components. Since 7i a , D a and E b , B b are real, the matrices 
Ql ba etc. are also real. The subject of our article will be the discussion of certain possible algebraic 
constraints between the components of the constitutive matrices Qi ba , < B ba ,€ b a ,D b a . The most no- 
torious one is the so-called Post constraint, which requires the sum of the traces of the matrices 
£ 6 a and D b a to vanish, i.e., <t a a + D a a = 0. There is an ongoing dispute mainly in the electrical 
engineering community, see Lakhtakia & Weiglhofer [12, 13, 14, 15, 16, 50, 51] and Lindell, Sihvola, 
Tretyakov, et al. [20, 21, 22, 23, 24, 38, 43, 49] whether this constraint is valid or not. We will 
address this question. 

The constitutive matrices are functions of space x a and time t in general. In other words, they 
are fields 2l 6a = 2l ba (x a ,t), etc. For applications to the case (i), to empty spacetime, they can be 
constants if a flat Minkowski spacetime and Cartesian coordinates are prescribed. However, if a 
Riemannian spacetime of general relativity is assumed or just a 4-dimensional differential mani- 
fold, then clearly the constitutive matrices are fields. Similarly, in case (ii), for material media, 
a homogeneous medium requires constant constitutive matrices, as soon as inhomogeneous media 
are allowed, the constitutive matrices are fields again. We may call them generalized permittiv- 
ity/impermeability fields. 

Before we turn to a 4-dimensional discussion of the Maxwell equations, we want to comment on 
the constitutive relations (5), (6). In the next section we will show that for reasons of 4-dimensional 
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covariance (5), (6) represent the appropriate constitutive relations. However, in the engineering as 
well as in the physics literature, we often find T> a and B a expressed in terms of and H^. If we 
exclude singular cases, i.e., 

det <B ^ , (7) 

then (5) can be multiplied with Q3 -1 , solved with respect to B a , and substituted into (6). In this 
way, we find, suppressing indices, 

V = (—21 + D*B _1 £) E + (2D*B _1 ) H , (8) 
B = (<B- 1 £)E+ (Q3" 1 ) H . (9) 

Thus, the the E'S-formalism in (5), (6) is transformed to the ETi-formalism of (8), (9). Both for- 
malisms are equivalent, provided a singular case can be excluded. The natural linear constitute 
relations are represented by (5), (6), as was shown by Post [34] and O'Dell [31], e.g., and as we will 
discuss in the next section. Hence we will concentrate on (5), (6) and will consider (8), (9) only as 
secondary. 



2 Electrodynamics in 4- dimensional spacetime 

The true structure of classical electrodynamics can be more clearly recognized if we present it in a 
4-dimensional way. If we define excitation, electric current, and field strength according to 



H = V - HAdt, J = p- j Adt , (10) 

F = B + E Adt, (11) 

then the 4-dimensional Maxwell equations (1) and (2) read 

dH = J, (12) 

dF = 0. (13) 



This shows explicitly that they are invariant under 4-dimensional coordinate/frame transformations 
and independent of the metric of spacetime, see [5]. If we require metric independence, we get a 
well-prescribed way of formulating the Maxwell equations. Since the metric of spacetime is the 
gravitational potential in general relativity, it is gratifying to know that there is a gravity-free way 
of formulating the Maxwell equations. Accordingly, the Maxwell equations (12), (13) are valid in a 
flat Minkowski spacetime (in any coordinates), in a curved Riemannian spacetime, and even in the 
Riemann-Cartan spacetime of the Poincare gauge theory of gravitation. Needless to say that this 
beautiful formalism has great practical value for applying electrodynamics in accelerated reference 
frames, for instance. 

The excitation H is a twisted 2-form, the field strength F an untwisted 2-form, see the discussion 
of Post in [36]. We decompose them in their components according to H = H^dx 1 A dx^ /2 and 
F = Fij dx l A dxi /2. Here x % are spacetime coordinates with = 0, 1, 2, 3. Analogously, the 

current J as a twisted 3-form decomposes as J = Jijk dx % A dx 1 A dx k /6. Thus, in the formalism of 
tensor analysis, (12) and (13) can alternatively be written as 

d k H lh = J\ (14) 
d {i F kl] = 0, (15) 



3 



with H lk := e tklm Hi m /2 and J 1 := e lklm Jkim/Q- Here e lklm = ±1,0 is the 4- dimensional contravari- 
ant Levi-Civita symbol. The covariant Levi-Civita symbol, which we will use below, is denoted by 
a circumflex: ej^ m — ±1,0. Incidentally, since no metric is available for the raising and lowering of 
indices at this stage, we have to distinguish these two symbols. 

In the constitutive law (5), (6), we assumed locality and linearity. In 4 dimensions this translates 
into 



H = k(F) , (16) 
where k is a local and linear operator. Accordingly, 

H H = g Ki i kl Fkl with Ki i H = ~ = ~ Ki i k ■ ( 17 ) 

Here Kif [x) is the twisted constitutive tensor of type [|]. It can characterize (i) either empty 
spacetime or (ii) a material medium. It has at most 36 independent components. In general, it is 
a field, that is, it can describe an inhomogeneous spacetime or an inhomogeneous medium. 

The relativistic invariance properties are crucial for deciding between the EB- and E^-formal- 
isms. From the point of view of a 4-dimensional spacetime, E and B form the invariant field 
strength F; likewise, T> and 7i form the invariant excitation H. Accordingly, the linear constitutive 
law (5), (6) is a relativistically invariant statement, as it is seen from (16). This law involves the 
constitutive matrices 21, 03, <£, and D that are combined into a covariant object - the constitutive 
tensor Ki kl . In contrast, the pairs T>, B and E, 7i do not form any 4-dimensional relativistic 
objects, and, as a result, the constitutive relation in the form of (8), (9) involves non-covariant 
objects with unclear physical meaning. Therefore, the EB-scheme of (5), (6) is, from a relativistic 
point of view, to be preferred against the £'?i-scheme of (8), (9), see O'Dell [31], Chap. 2, Sec. 1.2, 
pp. 23-24. If det ^ 0, both systems are formally equivalent; still, (8), (9), as a historical artifact, 
should be phased out from use. 

One can decompose the constitutive tensor Ki kl into its irreducible pieces under the linear group 
acting in the tangent and the cotangent space, respectively. Contraction is the only tool for such 
a decomposition. Following Post [35], we can define the contracted twisted tensor of type [\] 

:= K U kl , (18) 

with 16 independent components. The second contraction yields the twisted scalar function 

K ■= Kk k = K kl kl (19) 

(also called pseudo- or axial-scalar). The traceless piece 

t i h :=K i k -^K6t (20) 

has 15 independent components. These pieces (19) and (20) can now be subtracted from the original 
constitutive tensor. Then we find the following three irreducible pieces: 

K .M (l) K .kl , (2) kl , (3) kl 

= 0)^11 + 2^^ + 1^. (21) 
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By construction, ( > Kj-i j IS totally traceless: 

{1) K U kl = 0. (22) 

Thus, we split Ki kl according to 36 = 20 + 15 + 1, since the [|] tensor ^'k^ 1 is subject to the 16 
constraints (22) and hence carries 36 — 16 = 20 components. We call Kij kl the principal (or the 
metric-dilaton) part of the constitutive tensor. Without such a term, the propagation of electro- 
magnetic waves is ruled out, see the discussion subsequent to eq.(53). We further identify the two 
other irreducible parts ^Kij kl and ( 3 ^Kij kl as skewon and axion fields, respectively. Conventionally, 
the skewon and the axion fields are introduced with different factors by 

a '-=T 2 K - ( 23 ) 

If we substitute (21) into (17) i and use (23), then the general linear constitutive law can be 
written as 

Ha = \ {1) n t kl F kl + 2fi ]i k F j]k + a F tJ . (24) 

It has been shown recently [17] that taking the linear spacetime relation for granted, one can derive 
a Riemannian light cone provided one forbids birefringence in vacuum. However, if one allows for 
the skewon field, then the light cone is dissolved and we find birefringence and more general type 
of optical effects [30]. 

Let us now turn to some examples for case (i), the vacuum. 



Maxwell-Lorentz vacuum electrodynamics 

If we assume a Riemannian metric ^ with Lorentzian signature, then vacuum electrodynamics is 
recovered by 

(1) «y W = AoV=0 e ijmn g mk g nl , # = 0, a = 0. (25) 

Here Ao = \ATqJJm) is the vacuum admittance and g := det g k i- The principal part ^Ki kl corre- 
sponds to the Hodge star operator -k. The Maxwell equations can then be written as 1 

\ d*F = J, dF = 0. (26) 



Axion electrodynamics 

If we admit additionally the axion field a ^ 0, then H = X *F + ocF. This framework is called 
axion (Maxwell-Lorentz) electrodynamics, see Ni [25, 26, 27] and Wilczek [52], e.g.: 

\ d*F + (da) AF = J , dF = . (27) 

It is as if the current J picked up an additional piece depending on the gradient of the axion field. 
The real part of Kiehn's chiral vacuum [9] is a subcase, for a = const, of axion electrodynamics. 

Decomposing the inhomogeneous equation in time and space, we have (c is the speed of light 
and £ the 3-dimensional Hodge star) 

^-d^E = p-(da)AB, (28) 

\ cd^B -—d t ^E = j + (da)AE + (d t a)B. (29) 
c 

1 The inhomogeneous equation in indices reads Aq dj {yj—gg 3* Fki) = J 1 ■ 



5 



Possible experimental verifications were suggested by Sikivie and others, see [44, 45, 32]. 



Kij 



kl 



Pure axion field and the gyrator 

As a degenerate special case, we can also consider the pure ( "stand-alone" ) axion field with W 
(2) K .>* = o. Then, 

H = aF or H = -a E , V = aB, (30) 
and the Maxwell equations read 

{da)AF = J and dF = . (31) 

This is a special case of axion electrodynamics, namely (27) with Ao = 0. Here we have a case 
in which the constitutive matrix 03 in (5) vanishes. Consequently, det 03 = and the ET-C-version 
(8), (9) of the constitutive law loses its meaning. 

Seemingly the first person to discuss (and to reject) a constant pure axion field was Schrodinger 
[41], p. 25, penultimate paragraph, and, as a non-constant field, Dicke [3]. The framework (30), (31) 
corresponds to Tellegen's gyrator [46, 47] and to Lindell & Sihvola's perfect electromagnetic con- 
ductor (PEMC) [20, 21]. 

Lindell & Sihvola define their PEMC by using directly the relation H = aF. Tellegen [46, 47] 
considered a linear network of two terminal pairs (four-pole or two-port): 

vi = a n i 1 + a 12 i2, (32) 
v 2 = a 2 \ i\ + a 22 i2 ■ (33) 

The v's represent voltages, the i's currents, and the (frequency dependent) a's resistances. Since 

SI SI 

the field strengths are related to the voltage (E ~ V/m, B ~ Vs/m 2 ) and the excitations to the 

SI SI 

current (7i ~ A/m, V ~ As/m 2 ), our translation into the linear medium picture reads: 

E = a u V + a 12 H, (34) 
B = a 21 V + a 22 H. (35) 

Tellegen had a different translation, since he used the £7i-formalism of (8), (9). Tellegen [46], 
eq.(10), defined his gyrator by 

vi = —ai 2 , (36) 

v 2 = aii. (37) 

The gyrator "rotates" currents into voltages and vice versa. In our medium picture we have 

E = -aH, (38) 
B = aV, (39) 

i.e., the excitations are rotated into the field strengths. As a comparison with (30) shows, this 
corresponds to the axion field. Lindell & Sihvola [20] have shown that by using the .EH-formalism 
and by performing a suitable limiting transition (since det 03 = 0), one can also arrive at (38), (39). 
Consequently, it is established that the medium analogue of Tellegen's gyrator is the axion field - 
and both are equivalent to the PEMC of Lindell & Sihvola. 

The gyrator has been first realized for microwaves by the use of premagnetized ferrites. A 
corresponding figure with detailed explanations can be found in the review article of Hogan [6], 
Fig.12. 
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3 The constitutive tensor \ ljkl an d its properties 



We can raise the two lower indices of the constitutive tensor Ki kl with the help of the contravariant 
Levi-Civita symbol e^ mn . Thus, we can introduce the constitutive tensor density of spacetime, 

X l]kl --=\^ mn K mn k \ (40) 

which is equivalent to Kij kl but more widely in use. The notation x was chosen in order to conform 
with Post's convention [34]. The constitutive law can now be written in the conventional form as 

H ij = ^X ijkl Fki, (41) 

see Post [34], eq.(6.12). The constitutive tensor with all 36 components has also been discussed in 
an equivalent way by Lindell in his book [19], see also [42, 11, 43, 38, 23, 24, 20, 21]. 
We can decompose x^ kl irreducibly under the linear group, too. With the expressions 

( V W = \^ mn(X) ^mn k \ (42) 
{ 2 )yAi kl — 1 e ijrnn i 2 ) Kmn kl = — ^3 m [ k (43) 
(3) x ijkl = ijmn (3) ^kl = ijkl K ^ ^ 

2 X2 

we find 2 

^ijkl _ (l)^ijkl _j_ (2)^ijkl _|_ (3)^ijkl ^ ^) 
36 principal 20 skewon 15 axion 1 

The irreducible pieces carry the additional symmetries 

w x m = (i) x wj ^ (2) x ijM = _(2) x kuj ^ (3) x mi = (3) x [ijki] _ ^ 46 ) 

Clearly, the principal part with its 20 independent components is the one discussed by Post [34]. 
The skewon part with its 15 components vanishes if one assumes the existence of a Lagrangian 4- 
form from which the spacetime relation can be derived in full, see below. However, this assumption 
is not a necessary, as we have shown [5]. Finally, the axion piece with only 1 independent component 
is totally antisymmetric: 

W x ijkl = ae ijkl . (47) 

In Post's book [34], eq.(6.18), this piece is forbidden; therefore, Lakhtakia and Weiglhofer [14] called 
it the Post constraint: 

(3) x ijki _ q Qr a = (Post constraint) . (48) 



2 We want to underline that the irreducible decomposition of is algebraically similar the decomposition of a 
curvature tensor in a Ricmann-Cartan space (Ricmannian space with a metric compatible connection). Of course, 
only in a Ricmann-Cartan space the curvature tensor has 36 independent components and can carry a totally 
antisymmetric piece 

R [W], This piece, by the first Bianchi identity, is directly related to Cartan's torsion. But 
let us stress, these are just algebraic analogies, not more and not less. Prcmctric electrodynamics, as represented 
by eqs.(l) and (2), does neither couple to the metric nor to the connection of spacetime. Thus also torsion is not 
involved. Moreover, the spacetime relation eventually defines a metric, but that's all. 
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Post even carries on and requires c^a = 0, see his eq.(6.19). The reasons he gave for the latter 
constraint are incomprehensible to us. Therefore, we will not dwell on it any further. 

We want now to turn to criteria that will help us to decide which pieces of the constitutive tensor 
are reasonable from a physical point of view. Our kinematical analysis so far gave no convincing 
argument why any of its 36 components should not exist in nature. In the sequel, we will distinguish 
between case (i), the vacuum, and case (ii), material media. 

In electrodynamics, after formulating the Maxwell equations, one has to specify the energy- 
momentum density of the electromagnetic field. Also this can be done in a metric- free environment, 
see [5]. The energy-momentum 3-form in (holonomic) coordinates reads 

E i = ^{FAe i \H-HAe i \F) , (49) 

with Ci = di as the 4-dimensional vector basis of the tangent space. Let us denote the basis of 
the 4-dimensional 3-forms by e*. = €ki mn dx l A dx m A dx n /6. If we then decompose the energy- 
momentum 3-form with respect to e^, that is, Sj = % k e^, we find the conventional Minkowski 
energy-momentum tensor 

% k = ^6?F lrn H lm - F a H kl . (50) 

Of course, (49) and (50) are equivalent representations of the energy-momentum distribution. 

The representation (49) is most convenient for showing that the axion piece, even as field, doesn't 
carry electromagnetic energy-momentum. The constitutive law H = a(x)F yields F A ei\H = 
F A (ctF) = aF A &i\F. Together with the second piece in (49), this vanishes. For the principal 
piece it can be shown that it leads to positive energy, see [8] , also the skewon piece has non- vanishing 
contributions in general, see Table 1. We are talking here about the electromagnetic energy. If the 
axion field contributed an own kinetic term a la ~ g th (dia)(dk<y)/2 to the Lagrangian, then it would 
carry own energy. 

Table 1: Contribution of each irreducible part ( r 'x of the constitutive tensor density to the electro- 
magnetic energy-momentum, to the Lagrangian, and to the Tamm-Rubilar tensor density, see [5]. 
Our results are consistent with the findings of Kiehn [9, 10]. 



Irreducible 
part 


Contributes to 


energy-momentum 


Lagrangian 


TR tensor density 


principal ^'x 


yes 


yes 


yes 


skewon ^ 2 >x or 


yes 


no 


yes 


axion ^x or a 


no 


yes 


no 



This brings us to the Lagrangian of the electromagnetic field. The Maxwell equations (12), (13) 
can be derived from electric charge and magnetic flux conservation and the energy-momentum 
current (49) via the Lorentz force density. A Lagrangian is not needed in this procedure. The 
same it true for the local and linear constitutive law (in the sense of irreversible thermodynamics). 
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However, one can find an electromagnetic Lagrangian 4-form 

V=--HAF. (51) 

By substituting H = n(F) into the Lagrangian, the piece ^k(F) drops out because of the "anti- 
Onsager" symmetry (46)2- This shows that a possible skewon piece is related to dissipative effects. 
For material media this is certainly a possibility. For empty spacetime it may be a legitimate 
hypothesis. 

We will now turn to the propagation of the electromagnetic field (of "light"). In the geometrical 
optics limit, we find for the wave covector q the Fresnel equation [29, 5] 

G m {x)mMi = Q, (52) 
with the (metric-free) Tamm-Rubilar tensor density 

G ijkl (x) := j { e mnpq e rstu X mnr{l X jlps]k X l)qtu • (53) 

By simple algebra it is possible to prove that x^ 1 drops out from the Tamm-Rubilar tensor: 
gijki((i) x + (2) x + (3) x ) = gm^i) x + x (2)y Furthermore, G ijkl (^x) = 0. In other words, the 
axion part w^i*' = e % ^ kl a(x), in the geometrical optics limit, does not influence the propagation of 
light. For the skewon piece we find Q^ kl (^ X ) — 0- But even more far-reaching, by straightforward 
algebra, see [5], it can be shown that Q l ^ kl (^x + x) = 0- Consequently, for ^'x — 0, the 
Fresnel equation collapses and there is no orderly wave propagation in the geometrical optics limit. 
Therefore, the existence of conventional wave propagation requires ^x 0; that is, the principal 
part of the constitutive tensor must be non-vanishing. 

Although from Table 1 one may have an impression that the axion field is to a great extent 
irrelevant in experiments, it does affect physical observations. In particular, the important phe- 
nomenon of electromagnetic waves is certainly not exhausted by geometrical optics. Specifically, 
in the framework of axion electrodynamics (27) with a constant principal part Wx^ kl , Itin [7] (see 
also Carroll et al. [2]) has shown that the axion field does influence electromagnetic waves if one 
looks into exact solutions. Assuming the usual plane wave ansatz for the electromagnetic field 
strength F = f e lip , with the constant amplitude 2-form / and the phase function ip(x l ), we find a 
generalization of the Fresnel equation (52) in the form of [7] 

Q %3k \x) qWhQi ~ X <iW (diaO($a) W = • (54) 

Here qi = di<p is the wave covector, as usual. This extended Fresnel equation is no longer an 
homogeneous algebraic equation, and the axion field contributes explicitly to the second term. 
This necessarily leads, in the presence of the axion field, to the birefringence of the waves. 

As another manifestation of the axion field, let us consider the Lorentz force density j\ and the 
balance equation related to it (see [5], Sec. B5.1): 

fi = dEi + Xi, with X i :=-~(FAC i H-HA£ i F). (55) 

This is the force density that an electromagnetic field exerts on charges and currents. The axion 
field indeed drops out from the first term, which contains the electromagnetic energy-momentum 
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current (49). However, the second term is nontrivial. Using the elementary properties of the Lie 
derivative £j, we straightforwardly find for the constitutive law H = Xo*F + aF, 

Xi = -±(dia)FAF. (56) 

This additional force density arises from the axion field. 3 



4 Constitutive tensor split in space and time 

We will now come back to the constitutive relations (5), (6). We will again work with a 3-vector 
formalism. Again, we formulate it in a metric independent way. We decompose the electromagnetic 
field with respect to the 3-dimensional coframe $ a and the 2-form basis e& = e b cd$ c A d d /2: 

H = H a r, E = E a a , (57) 
V = V b e b , B = B b e b . (58) 

Then the constitutive relations (5), (6) can be put in the form of a matrix equation [5], 



V a I I yba 2J a J [ R b J \ / \ B 



« R T • (59) 



Moreover, by using the Levi-Civita symbol, see (40), we find for the x tensor density 

X) ^ ( B c t % ) ■ (60) 

If we compare (60) with (17)i, then, by straightforward algebra, the constitutive 3x3 matrices 
A, B, C, T> can be related to the 4-dimensional constitutive tensor density (40) by 

St 6 " = X° a °\ < B b a = \iac d hefX Cdef , (61) 

^6 = \i bcd X cd0 \ V a b = \e acd x" bcd - (62) 

Let us now employ the irreducible decomposition. For the principal part of k or x, respectively, 
we take the notation of Post [34], namely e ab = e ba (6 independent components), //~ 6 = [x^l (6 
components), and r ) a bl with 7 C C = (8 components). The skewon part we parametrize according to 



■ a " m (63) 
n b 



with, in 3 dimensions, the tensor skewon field s a b (9 independent components), the vector skewon 
m a (3 components), and the covector skewon n a (3 components). If we insert all of this in k, we 
eventually find 



l b a »a b l \ , ( s a b + 5 b a s c c -e abc m c \ , / 5 b a 



e ab j a b J V e n c s b a - 5^s c c J 1 V <5 b a 



« n xa ■ (64) 



principal part 20 comp. skewon part 15 comp. axion part 1 comp. 



3 There are proposals to combine the latter force term with the usual term tfiCj, thus arriving at a modified 
energy-momentum current [2, 7]. However, this only works under the additional assumption of the constancy of the 
gradient of the axion field, d{dia) = 0, and results in a non gauge invariant expression for the energy-momentum. 
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According to (59), we can evaluate the matrix elements: 

K a = (^-e abc m c )B b + {-^ a + s a b - 5 b a s c c ) E b - a E a , (65) 
V a = (e ab - e abc n e ) E b + ( 1 a b + s b a - 5 a b s c c ) B b + a B a . (66) 

Recall that e ab = e ba , n~ h l = fj,^ a , and 7 C C = 0. Incidentally, a is a 4-dimensional (axial) scalar, 
whereas s c c is only a 3-dimensional scalar. The cross-term 7% is related to the Fresnel-Fizeau 
effects. The skewon contributions m c ,n c are responsible for electric and magnetic Faraday effects, 
respectively, whereas the skewon terms s a b describe optical activity. 

5 Constitutive tensor and physics 

Up to here, we haven't said much about which piece of the constitutive tensor may be realized in 
nature. We just made a "kinematical" investigation of the most general local and linear constitutive 
tensor with its 36 independent components. From an analysis of the Tamm-Rubilar tensor we 
found that W x ijkl ^ is necessary for a non-degenerate propagation of electromagnetic waves in 
the geometrical optics limit. Only then the Tamm-Rubilar tensor doesn't vanish and the Fresnel 
equation exists. By the same method we find that, in the same limit, ( 3 > x l ^ kl ^ doesn't disturb 
the wave propagation at all. 

The skewon part ™x iJ ' fci , in its full generality, was discussed, amongst others, by Sihvola & 
Lindell [43, 19] and by us [5]. We studied already the possible effects of the skewon field on light 
propagation [30]. But we should mention Nieves and Pal [28] who introduced, besides the vacuum 
impedance 1/Ao and the velocity of light c, a third constant for the vacuum by means of a specially 
chosen skewon piece. Nieves and Pal assume (in our notation) 

s a b = i 5 b a , m a = , n a = (Nieves & Pal) . (67) 

This corresponds to a spatially isotropic skewon field. By (67), a constraint is given that restricts the 
class of allowed reference frames. Lakhtakia 4 calls it "isotropic chirality which is Lorentz reciprocal" . 
The Post constraint doesn't forbid it since the Nieves Sz Pal piece relates to the skewon part, whereas 
the violation of the Post constraint to the axion part. Lakhtakia points out that s 7^ is valid for 
certain materials, but not for vacuum. Furthermore, others (see [5]) put forward the hypothesis 
that there can be a fourth electromagnetic constant that makes the medium "nonreciprocal (but 
still isotropic)". It is, of course, the axion field. With (67), the constitutive law (65), (66) reads 

H a = - (s + a) E a + ^H a (E, B) , (68) 
V a = {i)x> a (E, B) - (s - a) B a , (69) 

where we used the notation ^TL := ^k(F) and ^T> := ^'k(F). We rewrite (68), (69) for the case 
of a Maxwell-Lorentzian principal piece in exterior calculus as 

H = [T 1 -B- (s + a)E } (70) 
V= e*-E -(s-a)B. (71) 

4 Private communication. 
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Note how the isotropic skewon field and the axion field act with different signs on E a and B a . 
Thereby they can be distinguished phenomenologically. We stress that the axion field a is a 
premetric, 4-dimensional scalar with twist and as such of fundamental potential importance. In 
contrast, the skewon field s is only a 3-dimensional scalar valid in a constrained class of reference 
frames. — 

What is then the the state of the Post constraint? Of course, the somewhat abstract concepts 
of a pure axion field, the Tellegen gyrator, and the PEMC of Lindell & Sihvola all violate the Post 
constraint, that is, a ^ 0, whereas the other pieces of the constitutive tensor vanish at the same 
time. 

One argument in favor of the Post constraint is the following (see [50]): With ^>H := ^k(F) 
and := ( 2 >n(F), the inhomogeneous Maxwell equation for the linear constitutive law reads 

dH = d( (1) H + i2) H) + (da)AF = J, (72) 

where we used already the homogeneous Maxwell equation dF = 0, see also (27). If the axion piece 
of the constitutive tensor is constant, i.e., da = 0, then a cannot couple to the Maxwell equations. 
However, this is only true for an infinitely extended material. As soon as the material has finite 
size, the axion piece a jumps at the surfaces of the material — and this jump of a couples to 
the Maxwell equations. In other words, the axion field can be seen along those jump surfaces by 
means of electromagnetic waves beyond the geometrical optics limit. In order to make this explicit, 
let us assume that a two-dimensional surface S divides the space into two parts. Let now both 
half-spaces be filled with homogeneous material media both of which are characterized by constant 
but unequal values of the axion field: a±, say, in the first material and «2 in the second one. Then, 
if we take as an example the pure gyrator case with the constitutive law (30), we find from the 
jump conditions on the boundary (see [5], Sec. B.4.3), 

S = -J°. (73) 

Here n a are the components of the covector normal to the surface S, and [a] := «2 — ct\ is the 
jump of the axion field on the boundary. Thus, even though the constant axion drops out from the 
Maxwell equations inside both homogeneous regions, it pops up in the form of the surface charge 
density p and the surface current density j a induced on the boundary S. Both are proportional 
to the jump of the axion [a]. In the more general case of axion electrodynamics, see (27), the 
left-hand sides in (73) will include the jump of the normal components of the electric excitation T> a 
and the jump of the tangential components of the magnetic excitation T~C a , respectively. Hence our 
argument can be generalized straightforwardly. 

This enables us to measure the value of the axion field of a homogeneous material medium. In 
this sense, the axion field of a material medium, our case (ii), always couples to the inhomogeneous 
Maxwell equation, and the argument mentioned above no longer applies. In the case of spacetime, 
also a curved one, the argument is correct if really everywhere da = 0. As soon as the axion field 
becomes space and/or time dependent, also in case (i) the quoted argument loses its meaning. 

For empty spacetime, our case (i), the Post constraint, according to experimental evidence, is 
fulfilled. However, the relaxation of the Post constraint leads to the concept of the axion field. 
According to our Table 1, this field is very elusive since it neither shows up in the electromagnetic 
energy-momentum tensor nor in the Tamm-Rubilar-tensor, i.e., it doesn't influence light propaga- 
tion in the limit of geometrical optics. Hence it is hard to detect it. A whole industry is at work to 
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find the axion field, see the new data on elementary particles [32], p. 389. So far without success. 
There is no theoretical reason known to us that would forbid the axion field. Quite the oppo- 
site. From a theoretical point of view, it appears as a quite natural extension of Maxwell-Lorentz 
electrodynamics. But, as mentioned, the axion hasn't been found so far. 

For material media, our case (ii), the situation is different. Gyrators have been engineered, 
see [6]. Consequently, there is little doubt that materials can be constructed that violate the 
Post constraint. The Tellegen model with parallel or antiparallel electric and magnetic dipoles 
looks reasonable, see Tretyakov et al. [49]. In the meantime, Tretyakov et al. [48] built an artificial 
Tellegen particle and verified its nonreciprocal magnetoelectric behavior. "Natural" material media 
exist that violate the Post constraint, namely C^Os and in Fe2Te Oq in static magnetic fields, see 
de Lange & Raab [18]. Accordingly, model calculations [33, 38] and experiments [18] (see also [37]) 
show that the Post constraint is valid for many materials; however, it is definitely violated in some 
rare cases. Thus, the Post constraint as a general dogma should be buried with all due honors. 
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